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Abstract. The package numericalsgps performs computations with and for numerical semigroups. 
Recently also affine semigroups are admitted as objects for calculations. This manuscript is a survey 
of what the package does^ and at the same time of the trending topics on numerical semigroups. 


1. Introduction 

The motivation for the implementation of numericalsgps was the lack of a package specific 
to make computations related to numerical semigroups. We had several functions implemented 
in distinct programming languages, each with its own interface, which made the communication 
between them very difficult. Thus we decided to unify all these procedures in a single package 
written in an appropriate language. The choice for the language was GAP ( |I29I| ), mainly because of 
the experience achieved by the first author. 

The first version of numericalsgps was released in 2005. Since then, the authors of this note 
have been adding new functionalities and replacing, when appropriate, algorithms with newer 
and faster ones. This makes the package to permanently reflect the state of the art in the area. 

Also some algorithms have different implementations that are used taking into account either 
the information stored in the object to deal with, or what packages has the user installed/loaded 
in his GAP installation. 

The current version is 1.0 and can be found in |26IJ ', the development version is available in 
https : / /bitbucket. org/gap-system/numericalsgps , The manual is over 90 pages long. 

The heart of this manuscript is Section [2j which consists of a brief description of the contents of 
the package. Aiming to make the paper self contained, we give definitions when necessary. The 
few examples given are simple illustrations that can guide the reader to produce his owns (assum¬ 
ing he has a working installation of GAP). We give pointers to the literature where one can find 
the implemented algorithms, which makes also the paper a kind of review of the computational 
procedures to deal with numerical and affine semigroups. 

The paper ends with a reference to interactions with other commutative algebra packages. The 
use of external software frequently improves very much the execution time of the functions offered 
for affine semigroups. 


2. Contents 

We briefly describe in this section the contents of the package, using as a guideline the chapters 
of the manual. 

2.1. Introduction. In the introduction of the manual, the basic definitions concerning numerical 
semigroups are given. The author interested in the topic can find all these definitions in [[40| . 
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2.2. Numerical semigroups. A numerical semigroup S is a submonoid of the set of nonnegative 
integers N under addition, and such that N \ S has finitely many elements. This section describes 
several ways to define a numerical semigroup. 

The elements in the set N \ S are usually called gaps, and its cardinality is the genus of S. We say 
that a gap g is a. fundamental gap of S if it is maximal in the set of gaps with respect to division, or in 
other words, g £ S, 2g e S and 3g e S. Gaps and fundamental gaps fully determine the semigroup 
S, and so, they can be used to describe a numerical semigroup. Functions are provided to deter¬ 
mine if a list of nonnegative integers is a list of gaps or fundamental gaps; and also procedures to 
define numerical semigroups by means of these lists. 

Since N \ S is finite, the maximum of Z \ S exists and it is known as the Frobenius number of S 
(there is actually a huge number of papers dealing with the computation/bounds of the Frobenius 
number of numerical semigroups; see for instance |38jF The conductor of S is just the Frobenius 
number of S plus one, and has the property that it is the least nonnegative integer c such that c -l- 
N c S. We call the elements in S less than or equal to the conductor the small elements of S. Clearly, 
the semigroup S is uniquely determined by its small elements. A procedure is implemented to 
check if a list of integers is the set of small elements of a numerical semigroup, and also a function 
to define a numerical semigroup if this is the case in terms of this list. 

If we take a closed interval I = [a/b,c/d] with a, b, c and d positive integers such that a/b < c/d, 
then the set IJjtgjvjjlNn/cf) is a numerical semigroup (and coincides with the set of all numerators of 
rational elements in I). It can be shown that this class of semigroups coincides with the numerical 
semigroups that consist of nonnegative integer solutions to inequalities of the form ax mod f < 
yx, which are known as proportionally modular numerical semigroups. If y = 1, then they are 
simply called modular. Hence we can also define a numerical semigroup in terms of the modular or 
proportionally modular inequality (giving a list with the parameters) or by an interval (providing 
its ends). Note that distinct intervals can yield the same numerical semigroup (and the same 
holds for proportionally modular inequalities). Membership to a numerical semigroup given by a 
proportionally modular inequality is trivial. Also specific fast algorithms exist for the computation 
of the Frobenius number if this is the case. For some kind of semigroups where testing being 
proportionally modular is fast, we perform this test and keep the inequality. 

Another way to uniquely determine a numerical semigroup is by any of its Apery sets of its 
nonzero elements. Let S be a numerical semigroup and let n 6 S \ {0}. The Apery set of n in S is the 
set {s€S|s-ngS}. This set has precisely n elements, one for each congruence class modulo n. 
Once we know an Apery set, the cost of membership to S is small, and also the Frobenius number 
and genus can be easily computed. Thus if the Apery set with respect to the least positive integer 
in S (its multiplicity) is computed, we store it as part of the object S. Many other invariants depend 
also on this specific Apery set as we will refer later. We provide a function to determine if a given 
list of integers is the Apery set of an element (the length of the list) in a numerical semigroup, and 
also to define a numerical semigroup by means of the Apery set. 

Apery sets and proportionally modular inequalities can be seen as particular cases of periodic 
subadditive functions. We say that / : IN —> IN is subadditive if /(O) = 0 and f{i + i) < f{i) + 
f(i) for all i,i e N. Associated to / we can define the semigroup of nonnegative integers x such 
that f{x) < X. This set is a numerical semigroup when / is periodic (with positive period). We 
represent a periodic function by the values on the integers less than the period, and thus by a list 
of nonnegative integers. We give a function to test if a list corresponds to a subadditive function, 
and if so, a numerical semigroup can be defined by using this list as argument. 

Let A be a nonempty subset of N. The monoid generated by A, denoted {A), is the set of all 
(finite) sums of elements of A. We say that A generates the numerical semigroup S if (A) = S. 
Observe that if this is the case, then the multiplicity of S must be in A, and whenever two gener¬ 
ators are congruent modulo the multiplicity, we do not need the largest one to generate the same 
semigroup. So we can always think of A to be finite (since all its elements can be chosen to be 
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incongruent modulo the multiplicity). Clearly, S is uniquely determined by any of its systems of 
generators. Among these, there is only one minimal with respect to inclusion (actually also with 
respect to cardinality) which is S* \ (S* + S*), where S* = S\ (0). The cardinality of this set is known 
as the embedding dimension of S. We give functions to define a numerical semigroup in terms of a 
generating set. 

gap> s:=NumericalSemigroup("interval",71/5,153/8); 

<Proportionally modular numerical semigroup satisfying 765x mod 10863 <= 197x > 
gap> t:=NumericalSemigroup(15, 16, 17, 18, 19, 29, 43); 

<Numerical semigroup with 7 generators> 

gap> s=t; 

true 

2.3. Basic operations. Among the basic operations of a numerical semigroup related to the con¬ 
tents of the preceding section, the package offers; computation of the multiplicity, generating 
system, minimal system of generators, small elements, gaps, embedding dimension. Apery sets, 
Frobenius number, conductor and fundamental gaps. 

Some functions have different methods depending on what is known about the semigroup. As 
an example, if the Apery set is known, the Frobenius number will be computed by using Selmer's 
formulas (see for instance [1401 Proposition 2.12]). 

Given S a numerical semigroup, we also give a procedure to list the first n integers in S, with n 
a positive integer. 

Associated to S we can define the partial order relation <5 on Z as follows. We write 

a<gbiib-aeS. 

The set of maximal elements of Z \ S with respect to this order are known as pseudo-Frobenius 
numbers (actually the Frobenius number is one of them), and their cardinality is the type of S. We 
provide functions to compute the pseudo-Frobenius numbers (that can be also obtained from the 
Apery sets) and the type of S. 

Let m be the multiplicity of S. Then the elements in the Apery set of m in S are Wj = kjm + i 
for i e {0,...,m- 1) with = 0 and [ki,...,km-\) fulfilling a set of inequalities (p2lj). In this 
way a numerical semigroup with multiplicity m corresponds with a point inside a polytope. We 
give a function that outputs the set of inequalities describing this polytope, and also to compute 
(k\,...,k„,-i), which are known as the Kunz coordinates of S. 

An element s e S is a minimal generator if S \ {sj is again a numerical semigroup. Hence the 
dual of this property could be an element g € S such that S U {g( is also a numerical semigroup. 
These elements are known as special gaps. We give a function to compute them, that can be used 
to compute oversemigroups of a given semigroup (Section [TSil . 

2.4. Presentations of a numerical semigroup. Let S be a numerical semigroup minimally gen¬ 
erated by {ni,...,Mg|. Then the monoid morphism cp : ^ S, (j 9 (ai,...,flg) = Yfi=i^ini is an epi- 

morphism, known as the factorization homomorphism of S. Consequently NVkertp is isomorphic 
to S, where keicp = {{a,b) e IN'^ x | (p{a) = (p{b)}. A presentation of S is a generating system of 
the congruence kercp. A minimal presentation of S is a minimal generating system of kercp (again, 
no matter if you think about minimal with respect to inclusion or to cardinality; both concepts 
coincide for numerical semigroups; see [1401 Chapter 7]). 

Minimal presentations can be computed from graphs associated to elements in the numerical 
semigroup. Let n be a nonzero element of S. We define the graph associated to n as the graph 
with vertices the generators n, e {Mi,...,ng[ such that n - n^ e S; and n^nj is an edge if n - (n, - 1 - 
nj) e S. There is a function to compute the graph associated to n. A minimal presentation is 
constructed from those graphs that are not connected (there are finitely many of them and can be 
found by using, once more the Apery set of the multiplicity). The elements having an associated 
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non connected graph are called Betti elements of S. A procedure to find the set of Betti elements of 
S is given in the package; and also to find a minimal presentation of S. 

Some numerical semigroups admit essentially a unique minimal presentation, in the sense that 
if a and r are two minimal presentations (and thus have the same cardinality), whenever (a, b) e a, 
either (a,b) 6 t or (b,a) e r (that is, unique up to permutation of the pairs of the presentation). 
In particular, generic numerical semigroups have unique minimal presentations ( [Il0| l. The semi¬ 
group S is generic if every pair {a,b) in a minimal presentation of S has the property that a-b 
has no zero coordinates. We give procedures to detect whether or not S is uniquely presented or 
generic. 

A straight generalization of the graph associated to n e S is the following: we can construct the 
simplicial complex of subsets A of {ni,...,Mg} such that n - ^ This set is known as the 

shaded set of n in S and has some nice properties associated to the generating function of S ([I45|l. 

The congruence ker (p is also a submonoid of xlN*^, which is generated by its nonzero minimal 
elements with respect to the usual partial ordering on x ( |I39[ Chapter 8]). If (a,b) is one 
of this minimal generators, then q}{a) = <p{b) 6 S is called a primitive element of S. These elements 
play an important role in factorization properties of S, and consequently we provide a function to 
compute them. 

gap> s:=NumericalSemigroup(5,7,9); 

<Numerical semigroup with 3 generators> 

gap> MinimalPresentationOfNumericalSemigroup(s); 

[ [ [ 0, 2, 0 ], [ 1, 0, 1 ] ], [ [ 4, 1, 0 ], [ 0, 0, 3 ] ], 

[ [ 5, 0, 0 ], [ 0, 1, 2 ] ] ] 

2.5. Constructing numerical semigroups from others. We have already seen that adding a spe¬ 
cial gap to a numerical semigroup produces a new numerical semigroup, and the same holds if 
we remove a minimal generator. The intersection of two numerical semigroups also produces a 
numerical semigroup. Functions performing these tasks are provided in numericalsgps. 

Let p be a positive integer and S be a numerical semigroup. The set S/p = {x 6 IN | px e S( is 
again a numerical semigroup, called the quotient of S by p. A function is given to compute this 
new semigroup. A kind of inverse is the notion of multiple of a numerical semigroup: given an 
integer a > 1 and a numerical semigroup S, then aS is a submonoid of N, but it is not a numerical 
semigroup. If we add to this set all the integers greater than or equal to a given positive integer, 
say b, then we obtain a numerical semigroup: aS U {&,—>}. If we start from N, and we repeat this 
operation several times, then we construct the set of what is known in the literature as inductive 
numerical semigroups (see for instance [|28| and the references therein) . 

For a numerical semigroup S the set of numerical semigroups T with S c T c IN is finite (the 
oversemigroups of S), since the genus of S is always finite by definition. We provide a function to 
compute the set of all oversemigroups of a given semigroup. Also there is a procedure to com¬ 
pute all numerical semigroups with given Frobenius number (this is done using the concept of 
fundamental gap as explained in [I43| |l and another function to compute the set of all numerical 
semigroups with given genus g (by constructing the tree of all numerical semigroups up to the 
level g). 

gap> s:=NumericalSemigroup(5,7,9); 

<Numerical semigroup with 3 generators> 

gap> Length(OverSemigroupsNumericalSemigroup(s)); 

15 

gap> Length(NumericalSemigroupsWithFrobeniusNumber(21)); 

1828 

We provide functions implementing the algorithms given in f27l| to compute the set of all nu¬ 
merical semigroups having a given set as set of pseudo-Frobenius numbers. 
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gap> pf := [13,24,25] ;; 

gap> NumericalSemigroupsWithPseudoFrobeniusNumbers(pf); 

[ ] 

gap> pf := [13,19,25];; 

gap> NumericalSemigroupsWithPseudoFrobeniusNumbers(pf); 

[ <Numericaf semigroup>, <Numerical semigroup>, <Numerical semigroup> ] 

2.6. Irreducible numerical semigroups. A numerical semigroup is irreducible if it cannot be ex¬ 
pressed as the intersection of two numerical semigroups properly containing it. This is equivalent 
to saying that it is maximal in the set of numerical semigroups with its same Frobenius number. 
Every numerical semigroup can be expressed (though not uniquely) as an intersection of irre¬ 
ducible numerical semigroups. We give a function to do this in numericafsgps (see [1401 Chapter 
3] for a description of the algorithm). 

We also give a procedure to compute all irreducible numerical semigroups with given Frobenius 
number: the procedure is based in [ llllj . This is actually equivalent to computing all irreducible 
numerical semigroups with given genus. This is due to the fact that if / is the Frobenius number 
of an irreducible numerical semigroup, then either g = {f + l)/2 or g = {f + 2)11, depending on the 
parity of /. 

A numerical semigroup S with Frobenius number / is symmetric if whenever X6 Z\S,/-X6 
S. The class of symmetric numerical semigroups coincides with that of irreducible numerical 
semigroups with odd Frobenius number. Irreducible numerical semigroups with even Frobenius 
number are called pseudo-symmetric. We give tests to detect if a numerical semigroup is in any of 
these classes. 

A particular class of irreducible numerical semigroups is the set of numerical semigroups with 
the least possible number of relations in its minimal presentations. These semigroups are called 
complete intersections, and it can be shown that every complete intersection numerical semigroup 
is either N or a gluing of two complete intersections (see for instance [1401 Chapter 8]). We say 
that S = Si -I- S 2 , with S a numerical semigroup and Si and S 2 submonoids of IN, is a gluing of Si 
and S 2 if gcd(Si)gcd(S 2 ) € Si n S 2 and gcd(Si) ^ I ^ gcd(S 2 ). We give procedures to detect if a 
numerical semigroup can be expressed as a gluing of two of its submonoids, and if it is a complete 
intersection. 

We also implement the procedures presented in ][3] to compute the set of all complete intersec¬ 
tion numerical semigroups with fixed Frobenius number (equivalently fixed genus, since we are 
still dealing with irreducible numerical semigroups). We present procedures to detect if a numer¬ 
ical semigroup is free (either IN or a gluing of a free numerical semigroup with a copy of N) and to 
calculate all free numerical semigroups with fixed Frobenius number. The same is done for tele¬ 
scopic numerical semigroups (these are free numerical semigroups where the gluing is performed 
in the same order given by the generators) and numerical semigroups associated to irreducible 
planar curve singularities (a particular case of telescopic numerical semigroups; see |3| for more 
details). 

A generalization of the concept of irreducible numerical semigroup is the following. We have 
seen that the genus g of an irreducible numerical semigroup S with Frobenius number / is either 
g = (f + l)/2ii f is odd (symmetric), or g = (/ - 1 - 2)12 if / is even (pseudo-symmetric). It turns out 
that the type of symmetric numerical semigroups is 1 and the type of pseudo-symmetric numerical 
semigroups is 2. So if S is an irreducible numerical semigroup with genus g, Frobenius number 
/ and type t, then g = (f -t- t)/2. We say that a numerical semigroup S is almost-symmetric if its 
genus is one half of its Frobenius number plus its type. We give a function to test if a numerical 
semigroup is almost-symmetric and include the procedure presented in 1141 1] to compute the set of 
almost symmetric numerical semigroups with fixed Frobenius number. 

gap> s:=NumericalSemigroup(5,7,9); 
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<Numerical semigroup with 3 generators> 
gap> Decomposelntolrreducibles(s); 

[ <Numerical semigroup>, <Numerical semigroup> ] 
gap> List(last,MinimalGeneratingSystem); 

[ [ 5, 7, 8, 9 ], [ 5, 1, 9, 11 ] ] 

gap> Length(TelescopicNumericalSemigroupsWithFrobeniusNumber(101)); 

86 

gap> Length(AlmostSymmetricNumericalSemigroupsWithFrobeniusNumber(31)); 

1827 

2.7. Ideals of numerical semigroups. A nonempty subset 7 of Z is a relative ideal of a numerical 
semigroup S if 7 + S c 7 and there exists d 6 Z such that d +1 c S (the concept of relative ideal 
corresponds to that of fractional ideal in domains). Every relative ideal 7 of S can be expressed in 
the form 7 = {ii,.. .,i„} + S for some integers ij. The set {li,. is a generating set of the ideal, and 
it is minimal if no proper subset generates the same ideal. 

gap> s:=NumericalSemigroup(3,4,5); 

<Proportionally modular numerical semigroup satisfying 5x mod 15 <= 2x > 
gap> 5+s; 

<Ideal of numerical semigroup> 
gap> [-1,2]+s; 

<Ideal of numerical semigroup> 
gap> MinimalGeneratingSystem(last); 

[ -1 ] 

We provide functions for computing the small elements of an ideal (the definition is analogous 
to that in numerical semigroups). Apery sets (and tables; see [I20||l, the ambient numerical semi¬ 
group, membership, and also some basic operations as addition, union, subtraction (I - J = {z e 
Z I z -I- / c 7)), set difference, multiplication by an integer, translation by an integer, intersection, 
blow-up (UneiN^^ “ ^rid >f-closure with respect to a family of ideals (|l44IJl. 

Numerical semigroups are “local" in the sense that there is a unique maximal ideal: the set 
of nonzero elements of the semigroup. Also there exists a canonical ideal, which for a numerical 
semigroup S with Frobenius number / is defined as{z€Z|/-zgS} (see for instance [[8|1. 

The Hilbert function associated to an ideal 7 of a numerical semigroup S is the function that 
maps every n 6 N to n7 \ (n -i- 1)7. The reduction number of 7 is the least positive integer n such 
that min(7) + nl = {n + 1)7. We give functions to compute the reduction number and the Hilbert 
function associated to an ideal. Also we give a procedure that computes the microinvariants of a 
numerical semigroup which are used to determine if the graded ring associated to the semigroup 
ring 7C[[S]] is Cohen-Macaulay (see HI). 

Finally we give a function to test if a numerical semigroup is a monomial semigroup ring fol¬ 
lowing |35l |. A numerical semigroup S is said to be monomial if for any ring R with K c Rc 77 [[xj 
and such that the algebraic closure of R is 7C[[x]] (77 a field with characteristic zero) and v(7?) = S, 
we have that 7? is a semigroup ring. Here, v denotes the usual valuation. 

2.8. Numerical semigroups with maximal embedding dimension. Recall that the embedding 
dimension of a numerical semigroup is the cardinality of its unique minimal generating system. 
Clearly, two minimal generators cannot be congruent modulo the multiplicity of the semigroup 
(the least positive integer in the semigroup). As a consequence, the embedding dimension is at 
most the multiplicity of the semigroup. Thus we say that a numerical semigroup S has maximal 
embedding dimension if its embedding dimension equals its multiplicity. 

The set of maximal embedding dimension numerical semigroups with fixed multiplicity, say m, 
is closed under intersection, and also if S {0} U (m-i-IN), then S U {/} is also a maximal embedding 
numerical semigroup, with / the Frobenius number of S. This in particular implies that if we are 
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given a numerical semigroup that is not of maximal embedding dimension, we can consider the set 
of all maximal embedding dimension numerical semigroups with its same multiplicity containing 
it, and then the intersection of all of them, obtaining in this way the maximal embedding dimen¬ 
sion closure of the given semigroup. Following this idea one can define the concept of minimal 
generators with respect to this class: the elements in the semigroup so that the closure of them 
yields the given semigroup. These elements are precisely the elements x in a maximal embedding 
dimension numerical semigroup S (together with the multiplicity) such that S \ jx) is a maximal 
embedding dimension numerical semigroup. 

gap> s:=NumericalSemigroup(3,5,7);; 

gap> MinimalMEDGeneratingSystemOfMEDNumericalSemigroup(s); 

[ 3, 5 ] 

We also give functions to compute the maximal embedding dimension closure of an arbitrary 
numerical semigroup. 

If S is a numerical semigroup with multiplicity m, then S has maximal embedding dimension 
if and only if for every x,y e S \ [0], x + y - m e S. A natural generalization of this pattern is 
the following. We say that a numerical semigroup S is Arf if for any x,y,z e S with x > y > z, 
then x + y - z e S. Clearly, every Arf numerical semigroup has maximal embedding dimension. 
Also, the class of Arf numerical semigroups is closed under finite intersections and the adjoin of 
the Frobenius number (of course if we are considering semigroups other than IN). Thus the class 
of Arf numerical semigroups is a Frobenius variety ( [1401 Chapter 6]). Again, it makes sense to 
talk about minimal generators with respect to this class, and also about the Arf closure of a given 
numerical semigroup (the intersection of all Arf numerical semigroups containing it). We give 
functions computing both things: Arf minimal generating sets and Arf closures. Also we provide 
a method to detect if a numerical semigroup is Arf, and a procedure that calculates the set of all 
Arf numerical semigroups with given Frobenius number. 

Finally, we consider in this section the class of saturated numerical semigroups, which turns 
out to be again a Frobenius variety (closed under intersections and the adjoint of the Frobenius 
number). A numerical semigroup is saturated if for every s,Si,...,s,- e S with s,- < s for all i and 
every zi,...,Zr e Z such that ZiSi -i- -i- z^s,- > 0 one gets s + ZiSi -i- -i- Zj-s^ € S. We provide for 

saturated semigroups the analogous functions that we described in the above paragraph for Arf 
semigroups. 

2.9. Nonunique invariants for factorizations in numerical semigroups. Let S be a numerical 
semigroup minimally generated by Recall that we defined a monoid epimorphism 

in Section [2~4l cp : IN'^ ^ S, (j9(ai,...,flg) = aiUi -t- -h Observe that for s e S, Z(s) = <p“^(s) 

collects the different expressions of s in terms of the generators of S. Thus we say that Z(s) is the 
set of factorizations of s. The cardinality of Z(s) is usually known as the denumerant of s. We use 
RestritcedPartitions to compute the set factorizations of s. 

The number of connected components of the graph associated to s € S (Section 12.411 coincides 
with the number of connected components of the graph with vertices given by Z(s) and zz' is an 
edge provided that z-z' ^0. 

Given z = {zi,...,zfi a factorization of s e S, we write |z| to denote the length of z, |z| = Zi + —i-Zg. 
The maximal denumerant of s is the number of factorizations of s with maximal length. Even 
though the denumerant is not bounded while s increases in S, the maximal denumerant is finite 
and can be effectively computed ( [ll4l| L We include this algorithm in the package as well as tests 
for supersymmetry and additiveness (see | jl4IJ for details). 

Let A be the Apery set of in S. A subset L of is an Lshape associated to S if (1) L c Z(A) 
(the set of factorizations of the elements in A), (2) for every ae A, #(L fi Z(a)) = 1, and (3) for every 
Z 6 L, if Z' 6 IN'^”^ is such that Z' < Z, then V e L. These sets give information on the factorizations on 
numerical semigroups (|[2|), and this is why we have included a procedure to compute them. 
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The set of lengths of factorizations of s is always finite (due to Dickson's lemma) and conse¬ 
quently we can write it as {/i < •• • < If]. The set {?2 - h>h k-i] is the Delta set associated 

to s. The Delta set of S is the union of all the Delta sets of s. This set is finite, and its maximum is 
achieved in one of the Betti elements of S ([I17|J). 

The elasticity of an element s is the ratio between the maximal and minimal lengths of fac¬ 
torizations of s. It was introduced to measure how far is a domain from being half-factorial (all 
factorizations of all the elements have the same length). No numerical semigroup other than IN 
is half-factorial, which is a unique factorization monoid. We give a procedure to compute this 
invariant. 

Given z = (zi,...,Ze),z' = (zp...,z') € Z(s), we denote by z Az' = (min(zi,z)),.. .,min(z£,,z')), which 
corresponds to the "common part" of these factorizations. The distance between z and z' is the 
d(z,z') = max(|z - z A z'|, |z' - z A z'|). The catenary degree of s is the least positive integer such that 
for any two factorizations of s, there exists a chain of factorizations such that the distance between 
two consecutive factorizations is bounded by this integer. The catenary degree of S is defined as 
the supremum of the catenary degrees of its elements. This supremum is reached in one of its 
Betty elements ([I16j[l. We give procedures to compute the catenary degree of a set of factorizations 
and of a numerical semigroup. Also other variants of catenary degrees are included: adjacent, 
homogeneous, equal or monotone catenary degree (see | j33[ I32j| l. For the homogenization of a 
numerical semigroup we offer a series of auxiliary functions. 

The tame degree of s € S is the least positive integer t such that for every factorization z of s and 
every integer i e {!,...,e} such that s-n^ e S, there exists another factorization z' of s with nonzero 
ith coordinate and such that the distance to z is less than or equal to t (there exists a factorization in 
which n,- is involved at a distance at most t). The tame degree of the semigroup S is the supremum 
of all the tame degrees of its elements, and it is reached in one of its primitive elements (also in an 
element with associated noncomplete graph). We give functions to compute the tame degree of a 
set of factorizations and that of the semigroup. 

Recall that associated to the numerical semigroup S, we can define the partial order on Z, a <g b 
if b - a e S. Thus (Z, < 5 ) is a poset, and one can define the Mobius function associated to it. We 
implement the procedure presented in [II 5| . 

The last invariant we give procedures to compute is the (U-primality, which determines how far 
is an element from being prime. The co-primality of s 6 S is the least positive integer co such that 
whenever s <§ JiaeA^^ with A c S finite, there exists Q c A with #Q < co such that s <5 JiaeO^- 
Clearly, if the omega primality is one, then the element is prime, if we look at <5 as a division.The 
ct>-primality of the semigroup is the maximum of the cu-primalities of its minimal generators. 
Initially we used the algorithm presented in | |10IJ . Now we use a faster procedure implemented by 
C. O'Neill (see Section [2TT3|. 

gap> 1:=FactorizationsIntegerWRTList(100,[10,11,13,15]) ; 

[ [ 10, 0, 0, 0 ], [ 1, 7, 1, 0 ], [ 3, 4, 2, 0 ], [ 5, 1, 3, 0 ], 

[ 0, 2, 6, 0 ], [ 3, 5, 0, 1 ], [ 5, 2, 1, 1 ], [ 0, 3, 4, 1 ], 

[ 2, 0, 5, 1 ], [ 7, 0, 0, 2 ], [ 0, 4, 2, 2 ], [ 2, 1, 3, 2 ], 

[ 0, 5, 0, 3 ], [ 2, 2, 1, 3 ], [ 4, 0, 0, 4 ], [ 1, 0, 0, 6 ] ] 

gap> TameDegreeOfSetOfFactorizations(1); 

5 

gap> CatenaryDegreeOfSetOfFactorizations( 1 ); 

3 


2.10. Polynomials, formal series and numerical semigroups. Let S be a numerical semigroup. 
The Hilbert series (not to be confused with the ITilbert function in Section [221 is the formal series 
Hs(^) = LseS Clearly = 1/(1-x) = Ls€iN\s Hence Ps(s) = 1 + (^-1 )Es6]N\s = 
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(1 -x)H 5 (x) is a polynomial, which we call the polynomial associated to S (see |I36I| ). We provide 
functions to compute both the polynomial and Hilbert series of a numerical semigroup. 

It turns out that when S is a complete intersection, the polynomial associated to S has all its 
roots in the unit circumference (and zero is not a root, which by Kronecker's lemma means that 
all the roots are in the unit circle, or equivalently, it is a product of cyclotomic polynomials). We 
give functions to determine if a monic polynomial with integer coefficients has all its roots in the 
unit circle, and to do this we need two auxiliary implementations: that of being cyclotomic and 
the computation of the Graeffe polynomial (see [IlSIJ for details). A numerical semigroup is said to 
be cyclotomic if its associated polynomial has all its roots in the unit circle. 

Symmetry (see Section can also be characterized in terms of the associated polynomial: a 
numerical semigroup is symmetric if and only if its associated polynomial is self-reciprocal (a 
palindrome if we look at the coefficients). 

Let K be an algebraically closed field. And let / e K[x,y] represent an irreducible curve with 
one place at infinity. Take g e K[x,y] and set int(/,g) = dimf^{K[x,y])/(f,g). Then the set {int(/,g) | 
g € (/)( is a numerical semigroup. We give a procedure to implement it (see jll). This kind of 
semigroups are generated by what is called a 6 -sequence. There is a function to compute all 6- 
sequences with fixed Frobenius number (equivalently genus since these semigroups are complete 
intersections and thus symmetric). Also associated to any 6 -sequence there is a "canonical" planar 
curve, and we offer a method to compute it. 

Let F be a set of polynomials. Then the set of values (respectively degrees) of the series (respec¬ 
tively polynomials) in the algebra fC|[f]] (respectively J^[F]) is a submonoid of N. Under certain 
conditions it is a numerical semigroup, and we provide functions to compute it. Also to deter¬ 
mine a basis of the algebra i<C[[f ]] (or K[F]) such that the values (or degrees) minimally generate the 
semigroup of values of this algebra (see |]5|). 

gap> t:=Indeterminate(Rationals,"t");; 
gap> 1:=[t"4,t"6+t"7,t"13]; 

[ t^4, t^7+t^6, tH3 ] 

gap> SemigroupOfValuesOfCurve_Local(1); 

<Numerical semigroup with 4 generators> 
gap> MinimalGeneratingSystem(last); 

[ 4, 6, 13, 15 ] 

gap> SemigroupOfValuesOfCurve_Local(l,"basis"); 

[ t'4, t^7+t^6, tH3, tH5 ] 

2.11. Affine semigroups. An ajfine semigroup is a finitely generated submonoid of N” for some 
positive integer n. In the package, affine semigroups can be defined by means of generators, as the 
set of elements in the positive orthant of a subgroup of Z” (full semigroups) or as the set of ele¬ 
ments in the positive orthant of a cone (normal semigroups). Our intention is to provide as many 
functions as possible for affine semigroups as we offer for numerical semigroups. Along this line, 
we present methods for membership, computing minimal presentations, determine gluings, Betti 
and primitive elements, and the whole series of procedures for nonunique factorization invariants 
(an overview of the existing methods for the calculation of these invariants can be found in [I30IJ L 
New procedures are now under development based in Hilbert functions and binomial ideals ( [I37I] L 

As an example, let us do some computations with G n N^, where G is the subgroup of with 
defining equations x+y = 0 mod 2 and x+z = 0 mod 2 (this is actually the block monoid associated 
to Z^; see |33l | for the definition of block monoid). 

gap> a:=AffineSemigroup("equations",[[[1,1,0],[0,1,1]],[2,2]]); 

<Affine seniigroup> 

gap> GeneratorsOfAffineSemigroup(a); 

[ [ 0 , 0 , 2 ], [ 0 , 2 , 0 ], [ 2 , 0 , 0 ], [ 1 , 1 , 1 ] ] 
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gap> OmegaPrimalityOfAffineSemigroup(a); 

3 

gap> BettiElementsOfAffineSemigroup(a); 

[ [ 2 , 2 , 2 ] ] 

2.12. Random. Based on the the methods provided hy GAP to create "random" objects, we pro¬ 
vide some functions for "random" affine numerical semigroups. These are particularly useful to 
produce examples. Furthermore, they are extensively used each time new algorithms are imple¬ 
mented and tests need to be made. 

gap> l:=List([1..20], _->RandomNumericalSemigroup(5,200));; 

gap> Is:=Filtered(1, s-> 1+FrobeniusNumber(s)=Genus0fNumericalSemigroup(s)*2);; 
gap> Fist(ls,MinimalGeneratingSystem); 

[ [8, 103 ], [ 25, 109 ], [ 35, 57, 125 ], [ 3, 52 ], [ 15, 170, 178 ], 

[ 3, 145 ], [ 21, 68, 153 ] ] 

2.13. Contributions. There is a special section devoted to contributions. So far we are happy 
to count with functions implemented by A. Sammartano and C. O'Neill (apart from those co¬ 
implementations with J. I. Garcia-Garcia and A. Sanchez-R.-Navarro). 

The functions implemented by Sammartano are mainly focused on deriving properties of the 
semigroup algebra k[[S]] and its associated graded algebra from properties of the numerical semi¬ 
group S. He offers procedures to determine purity and M-purity of S ([I13jl. Buchbaum ([I21]l. 
Gorenstein ([1221]) and complete intersection (|]24|) property for the graded algebra; some special 
shapes of the Apery sets (a, /I and y-rectangular, see [I23||h and the type sequence of a numerical 
semigroup f|[7|1. 

O'Neill on his side offers methods dealing with non unique factorization invariants: factoriza¬ 
tions, (u-primality and Delta sets for a list of elements in a numerical semigroup. Delta sets for the 
whole semigroup, and periodicity for the Delta sets ( f9 |). 

3. Interaction with other packages 

Since the first release of the package many other packages have come into scene (some still under 
development). We have tried to take advantage of this. Dealing with affine and numerical semi¬ 
groups translates in many cases to computing nonnegative integer solutions of linear Diophantine 
equations or Grobner basis calculations of binomial ideals. Hence the interaction with singular 
(|25ljl, Normaliz ( [I12IJ 1 and 4ti2 ([|T|) was a step forward for us. For singular there are several op¬ 
tions to consider: [1191] , | f6| and Singularinteface https ://github. com/gap-system/Singular Interface 
As for 4ti2, we can use [I34| and [13If . which is under development. Finally there is an interface for 
Normaliz that can be found in https ://github. com/f ingolf in/Normalizinterf ace. We have 
implemented different methods for each procedure depending on which of the above packages 
the user has loaded/installed. 
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